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We apply the integrated perturbation theory [ 1 ] to evaluate the scale-dependent bias in the presence of pri- 
mordial non-Gaussianity. The integrated perturbation theory is a general framework of nonlinear perturbation 
theory, in which a broad class of bias models can be incorporated into perturbative evaluations of biased power 
spectrum and higher-order polyspectra. Approximations such as the high-peak limit or the peak-background 
split are not necessary to derive the scale-dependent bias in this framework. Applying the halo approach, previ- 
ously known formulas are re-derived as limiting cases of a general formula in this work, and it is implied that 
modifications should be made in general situations. Effects of redshift-space distortions are straightforwardly 
incorporated. It is found that the slope of the scale-dependent bias on large scales is determined only by the 
behavior of primordial bispectrum in the squeezed limit, and is not sensitive to bias models in general. It is 
the amplitude of scale-dependent bias that is sensitive to the bias models. The effects of redshift-space dis- 
tortions turn out to be quite small for the monopole component of the power spectrum, while the quadrupole 
component is proportional to the monopole component on large scales, and thus also sensitive to the primordial 
non-Gaussianity. 

PACS numbers: 98.80.-k, 98.65.-r, 98.80.Cq, 98.80.Es 



I. INTRODUCTION 

The primordial non-Gaussianity is a useful indicator in 
searching for the generation mechanism of density fluctua- 
tions in the universe. While the primordial non-Gaussianity 
is small in the simplest model of single-field slow-roll in- 
flation, there are various other inflationary models which 
predict fairly large primordial non-Gaussianities (for review, 
see Refs. [2, 3]). Alternative models without inflation can 
also produce large non-Gaussianities (see, e.g., Ref. [0]). 
Accordingly, constraining or detecting the primordial non- 
Gaussianity has a substantial importance in studying very 
early stages of the universe. 

The large-scale structure (LSS) of the universe has been 
one of the most important ways of constraining cosmological 
models in general. In recent years, it was found that the pri- 
mordial non-Gaussianity produces the scale-dependent bias 
in the LSS JD-0]. Although the scale-dependent bias from 
the primordial non-Gaussianity mainly appears on very large 
scales, the form of the scale dependence does not receive gen- 
eral relativistic corrections even on scales larger than the Hub- 
ble radius n. 

The scale-dependent bias as a method to constrain the pri- 
mordial non-Gaussianity has already been applied to observa- 
tional data of galaxies and quasars 1^41211. The constraints de- 
rived from the scale-dependent bias are competitive with the 
measurements in the cosmic microwave background (CMB). 
A hint of a positive value of local-type non-Gaussianity pa- 
rameter/NL was indicated by analyses of the scale-dependent 
bias II 1211 . which could have profound implications for mod- 
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els of the early universe. However, more detailed analyses 
with large galaxy surveys are necessary to derive conclusive 
results. 

On the theoretical side, analytic expressions for the scale- 
dependent bias has been only approximately derived. So far, 
at least three kinds of derivation are known. One is based 
on the method of peak-background split [l3l - fl8ll . and 
the other is based on the statistics of high-threshold regions 
la, 1 171 1 19l 12011 . It is also shown that the scale-dependent bias 
gen erally appears in phenomenological models of local bias 
ll2ll l22ll . Because those theoretical derivations are approx- 
imate and not exact, t hey should have been compared with 
numerical simulations lu3Tl27ll . The scale-dependent bias is 
qualitatively confirmed by simulations, although the detailed 
amplitude of the analytic predictions needs to be modified 

ei 11 a 

The purpose of this paper is to provide more precise and 
more general analytic expressions for the scale-dependent bias 
in the presence of primordial non-Gaussianity. Evolutions of 
density fluctuations on sufficiently large scales are expected to 
be described by the nonlinear perturbation theory. However, a 
consistent inclusion of the general form of bias in the nonlin- 
ear perturbation theory had not been clear until recently. 

The integrated perturbation theory (iPT) [Ql is the formal- 
ism in which a broad class of biasing scheme can be con- 
sistently included on a general ground. The local Eulerian 
biasing scheme has been frequently adopted in attempts of 
incorporating the bias into the nonlinear perturbation theory. 
However, the local Eulerian bias is not consistent with the 
nonlinear dynamics in general 01 HH, because the nonlinear 
evolutions are nonlocal phenomena. In the formalism of iPT, 
generally nonlocal biasing either in Eulerian and Lagrangian 
spaces can be consistently included in the nonlinear pertur- 
bation theory. The effects of primordial non-Gaussianity and 



2 



redshift-space distortions are naturally incorporated into the 
formalism. 

Consequently, it is a straightforward application of iPT to 
making predictions of scale-dependent bias in the presence of 
primordial non-Gaussianity. In this paper, we first present the 
most general prediction of iPT for the scale-dependent bias, 
which can be applied to almost any type of primordial non- 
Gaussianity and to almost any model of bias, as long as we 
consider the regime where the perturbation theory applies. We 
then consider popular models of primordial non-Gaussianity, 
i.e., local-, equilateral-, folded-, and orthogonal-type non- 
Gaussianities. Asymptotic behaviors of scale-dependent bias 
on large scales, which have been derived in limited cases in 
the literature, are re-derived in general situations without re- 
sorting to specific forms of bias. When the halo model of bias 
is adopted, the previously known formulas of scale-dependent 
bias are re-derived by taking appropriate limits of our gen- 
eral formula. In course of derivation, nonlocality of bias turns 
out to be important. We also show that previous formulas de- 
rived from the peak-background split are only consistent with 
the Press-Schechter mass function. When the mass function 
deviates from the Press-Schechter form, our general formula 
predicts that those previous formulas of scale-dependent bias 
should be corrected. 

The main purpose of this paper is to show how the iPT 
can be applied to making theoretical predictions of the scale- 
dependent bias, and to give lowest-order calculations with pri- 
mordial bispectra, when the bias is given by a simple, nonlocal 
model of halo bias. With the iPT formalism, predicting the 
scale-dependent bias is straightforward once the bias model 
is given. Theoretical uncertainties in predicting the scale- 
dependent bias are reduced to those of the bias model. Ac- 
curate modeling of biasing is actively studied in recent years. 
Once we have an accurate model of galaxy bias which is 
generally nonlocal, the iPT immediately gives predictions of 
scale-dependent bias with least approximations. 

This paper is organized as follows. In Sections HIUlVI an- 
alytic derivations of scale-dependent bias in real space are 
presented in order. The general formula of the biased power 
spectrum in real space with primordial bispectrum is derived 
by the lowest-order iPT in Sec. [TT] Large-scale limits of the 
scale-dependent bias in concrete models of primordial non- 
Gaussianity are generally investigated in Sec. |III]without as- 
suming the models of bias. In Sec.[lV] we consider the shapes 
of renormalized bias functions which are needed in iPT. We 
generalize the previous results to include the effects of the 
smoothing function in the halo model of bias. Numerical eval- 
uations of the analytic formula are presented and compared 
with previous predictions in Sec. [V] In Sec. [VI] we general- 
ize our formula to include the redshift-space distortions. In 
Sec. lVIII we summarize our results. 

In plotting the figures of this paper, we adopt a cosmolog- 
ical model of flat curvature with parameters £? m o = 0.275, 
Q\o = 1 — Q mQ = 0.725, 0,o = Ono - Ao = 0.046, n s = 0.96, 
cr g = 0.8, Ho = 70km/s/Mpc, where Q m o is the matter den- 
sity parameter, Ovo is the cosmological constant parameter, 
£2bo is the baryon density parameter, Q c q is the cold dark mat- 
ter density parameter, n s and cr 8 are respectively the spectral 



FIG. 1: The diagrammatic representation of the power spectrum in 
terms of multipoint propagators. Details of diagrammatic rules are 
described in Ref. Q||. 



index and the amplitude of primordial density fluctuations, 
and Ho is the Hubble's constant. 



II. INTEGRATED PERTURBATION THEORY AND 
SCALE-DEPENDENT BIAS FROM PRIMORDIAL 
NON-GAUSSIANITY 

We apply the formalism of iPT to derive a general for- 
mula of the scale-dependent bias from the primordial non- 
Gaussianity, without assuming the peak-background split or 
the high-peak limit. In the derivation, it is convenient to intro- 
duce the multipoint propagators of biased objects. The mul- 
tipoint propagators of mass and velocity fields are recently 
introduced in Ref. |29[ l30ll . and those of biased objects are 
introduced in Ref. The multipoint propagators of bi- 
ased objects X are defined by ensemble averages of functional 
derivatives 

/ « \ = (2n f-3n § 3 (k ^ n _ k}r <n) {k k) 

\66 L (k l )---66dk„)l D X ^ 

(1) 

where 5iS.k) is the Fourier transform of the linear density field, 
and 6x(k) is the Fourier transform of the number density field 
of the biased objects X in Eulerian space. On the right-hand 
side, k\...„ = ki+- ■ -+k„ and 5^(k) is the Dirac's delta function 
in three dimensions. 

In terms of the multipoint propagators, the biased power 
spectrum of objects X is represented by 

Px(k) = [r£\k)] 2 P h (k) 

+ 4 1) (*) j ^rf(k',k-k')B L (k,k',\k-k'\) 

+ •■■, (2) 

where is the linear power spectrum and Z?l(^i , ^2, £3) is 
the linear bispectrum. The diagrammatic representation of the 
above equation is shown in Fig.Q] The details of the diagram- 
matic rules are described in Ref. [ 1 ]. The multipoint propaga- 
tors are represented by grey circles and the linear poly spectra 
are represented by crossed circles. The multipoint propaga- 
tors include all the loop corrections at the corresponding ver- 
tices. For Gaussian initial conditions, the linear bispectrum is 
absent. In this paper, we consider only lowest-order contri- 
butions of the primordial non-Gaussianity which are linearly 
dependent on the bispectrum. Contributions from the trispec- 
trum and higher-order poly spectra, and terms like a product of 
power spectrum and bispectrum are all neglected. 

The multipoint propagators and llj? also contain con- 
tributions from primordial non-Gaussianity in general. How- 
ever, such contributions in Eq. (fj) are higher orders in the 
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FIG. 2: Diagrammatic representations of the first two multipoint 
propagators with renormalized Lagrangian bias functions in the 
lowest-order approximation. Details of diagrammatic rules are de- 
scribed in Ref. |l||. 



above sense. Accordingly, we do not have to consider the non- 
Gaussian corrections to evaluate the multipoint propagators of 
Eq. (Q~|i in the lowest-order approximation of this paper. 

In the iPT, a concept of renormalizing Lagrangian bias 
functions is introduced. The renormalized bias functions in 
Lagrangian space c\ is defined by 



J (27T) 



f \dd L m---d6 L (k n ) 



, (3) 



where 8^(k) is the Fourier transform of the number density 
field of biased objects X in Lagrangian space. An equivalent 
way of defining the renormalized bias functions is given by 



8 n 8\(k) 



\6S L (k l )---S6 L (k„) 



(2n) 3 - 3n S 3 D (k h .. n - k)c^k } 



,*„). 
(4) 

In this form of definition, the similarity of the renormalized 
bias functions in Lagrangian space and the multipoint prop- 
agators of Eulerian space is obvious. The renormalized bias 
functions c\ can be seen as multipoint propagators of biasing 
in Lagrangian space. When the bias is not present, the num- 
ber density field of objects is uniform in Lagrangian space, 
#x = 0, and all the renormalized bias functions vanish. 

Applying the diagrammatic methods of iPT [1] in the 
lowest-order approximation, the first and second multipoint 
propagators in terms of renormalized Lagrangian bias func- 
tions are given by 



r^(k) = k-Li(k) + c\(k) 
r%\k u k 2 ) = k ■ L 2 (k u k 2 ) + Ik ■ Liiki)] [k ■ L 1 (k 2 )] 



(5) 



+ k ■ Li(ki)c\(k 2 ) + k ■ L x (k 2 )c\{k x ) 

+ cfa u k 2 ), (6) 

where L„ is the nth order perturbation kernel of displacement 
field, and k = k\ + k 2 in Eq. (0. The diagrams for the above 
expressions are shown in Fig. [2] The correspondences be- 
tween diagrams and equations are detailed in Ref. We 
omit partial resummations of the displacement field *F. If 
we leave these partial resummations, the right-hand sides of 
Eqs. © and © are multiplied by a factor U(k) = (e~ li v ). 
This factor is almost unity on sufficiently large scales. More- 
over, this factor is cancelled when we calculate the scale- 
dependent bias Ab(k) of Eq. d22b below. Therefore, we do 



not include this factor from the first place. One should note 
that Eqs. (0) and (0 are lowest-order results. For higher-order 
calculations, loop corrections to the tree-level expressions of 
Eqs. © and (O for multipoint propagators should be included 
as depicted in Fig. 17 of Ref. yj]. 

The perturbation kernels L„ are calculated by the La- 
grangian perturbation theory ll3l1 - l37l . In real space, the first 
two kernels are given by 



3 k 12 



L 2 (k u k 2 ) = -—j 



1 - 



k\ ■ k 2 

k\k 2 



and Eqs. and (O reduce to 

rf(k) = i + c\(k) 



r%\k 1 ,k 2 ) = F 2 (kuk 2 ) + \l + 



k\ ■ k-> 



c\(k 2 ) 



tl-h\c\(k l ) + c\(k l .k:) 

k 2 



where 



F2ikl kz) = B + (h + b] hih + 1 IhJ^ 

7 \k 2 k\ ) k x k 2 7 I k\k 2 , 



(7) 
(8) 

(9) 
(10) 

(11) 



is the second-order kernel of the Eulerian perturbation theory 
[HH. In the literature, an extra factor of 1/2! is frequently 
multiplied to the right-hand side of Eq. (fTTt . We find it is 
convenient not to include the factor for symmetrical reasons 
in diagrammatic methods 01 . 

In the lowest-order approximation, non-Gaussian contribu- 
tions to the renormalized bias functions can be neglected as 
described earlier. Consequently, we can consider that the first 
term of right-hand side of Eq. (O corresponds to the Gaus- 
sian part P^(k) and the second term corresponds to the non- 
Gaussian part AP x (k): 

P°(k) = [li\k)] 2 P h (k), 



(12) 



AP x (k) = lf(k) f ^^(2n?8 D {k x + k 2 -k) 
J (2tt) j {2ny 



xr^(k u k 2 )B L (k,k u k 2 ). (13) 



One should note that the above distinction holds only in 
the lowest-order approximation. In general, there are non- 
Gaussian corrections to the first multipoint propagator T^ 1 ' 
beyond the lowest-order approximation. Those corrections 
in Eq. ( fT2l should be included in APx{k) if we consider the 
higher-order approximations. For example, the lowest-order 
non-Gaussian correction term is proportional to the primor- 
dial bispectrum. However, such correction terms are multi- 
plied by the power spectrum in Eq. (TTZV Accordingly, the 
orders of non-Gaussian corrections to the multipoint propa- 
gator in Eq. (fT2l are higher than those of Eq. (Q~3}, and can 
be neglected in the lowest-order approximation. For the same 
reason, we do not consider loop corrections to the Gaussian 
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part Py (k) which are extensively studied in previous papers 
ll36l l39M4lll with Gaussian initial conditions. Evaluations of 
the leading non-Gaussian part AP x {k) is the main subject in 
this paper. 

Substituting Eqs. (O and ( TTOt into Eqs. ( fT2l and ( fT3l > we 

have 



P°(k) = [bdk)} 2 Pdk), 
AP x (k) = bdk) [Qo(k) + Qdk) + Q 2 (k)] . 



(14) 
(15) 



where 



bi(k) 



1 + c\(k) 



corresponds to the Eulerian linear bias factor, and the func- 
tions Q„{k) are defined by 



Qn(k) 



f d 3 k' - 



(k,k')B L (k,k',\k-k'\), (17) 



where 



&>(*,*') = 



&(*,*') 



2k k' 4 k z k' z - (k ■ k') 
— 

2k- k' 



7 k' 2 \k-k'\ 2 



c\(\k-k'\), 



Q 2 (k, k') = c\{k',k-k'). 



(18) 

(19) 
(20) 



When the bias is local, the linear bias factor b\{k) is scale- 
independent. In general, the scale-dependent bias emerges 
from nonlocality even when the primordial non-Gaussianity 
is absent. 

When the scale-dependent bias Ab(k) from the primordial 
non-Gaussianity is defined by 



p x (k) = [bm + Ab(k)f pjk), 

and higher orders of Ab are neglected, we have 



Ab(k) = IbiQO 



( AP x (k) AP m (k) 
P%(k) P?Jk) ) 



(21) 



(22) 



where P^(k) is the power spectrum of mass without primor- 
dial non-Gaussianity and AP m {k) is the contribution from the 
primordial non-Gaussianity, so that P m (k) = P^(k) + AP m (k) 
is the total power spectrum of mass. These are given by sub- 

c\ = (b x 
APx(k), respectively: 



stituting c\ — = (b\ = 1) in the expressions of P^(k) and 



F%k) = P L (k), 
4Pm(k) = Qo(k). 



(23) 
(24) 



Substituting Eqs. {14}, (Q3), and (O into Eq. d22l) . we 
have 



Ab(k) 



1 



2P L « 



[(M*) + &(*)-<$•(*)&>(*)] . (25) 



The above expression for the scale-dependent bias is con- 
sistent to that derived from the cross power spectrum Pmx(k) 



between mass and objects X as shown below. In fact, the cross 
power spectrum is given by 



-<i) 



Pmxik) = r^(k)r^(k)P h (k) 
1 

+ - 

2 

1 

2 



f ^rf(k',k-k')B h (k,k',\k~k'\) 
J (2tt) j 

f ^rf(k',k-k')B L (k,k',\k-k'\) 

(26) 



where T^' is the nth multipoint propagator of mass. In the 
( ) lowest-order approximation, the first two propagators are ob- 



tained by substituting c\ — = in Eqs. (0 and ( fXOt : 

lt\k) = 1, 
I%\k lt ki) = F 2 (k u k 2 ). 



(27) 
(28) 



Substituting Eqs. ©, CGI, d27]) and ([28j into Eq. (|26]l, we 

have 



AP mX {k) = 1+b 2 m Qo(k) + + X -Q 2 (k). 

When the scale-dependent bias Ab(k) is defined by 

P m x(k) = [b 1 (k)+Ab(k)]P m (k), 
and higher orders of Ab are neglected, we have 
(AP mX (k) AP m (ky 



Ab(k) = fei(jfc) 



(29) 
(30) 

(31) 
(32) 



Substituting Eqs. (|23), {24}, (|29]) and {30]) into Eq. {32), we 
again have the same equation as Eq. (f25). 

Our result, Eq. (|25t , is a general formula and any approxi- 
mation other than the perturbation theory is not employed. In 
the literature, analytic formulas of scale-dependent bias from 
the primordial non-Gaussianity have been derived either in the 
high-peak limit or in the approximation of peak-background 
split. These approximate results are subclasses of our general 
formula, as we explicitly show in the rest of this paper. 



III. MODELS OF PRIMORDIAL NON-GAUSSIANITY 

A. Specific models of primordial bispectra 

Models of primordial non-Gaussianity are quite com- 
monly characterized by the primordial bispectra B& of gauge- 
invariant Newtonian potential <2> at the matter-dominated 
epoch. The linear density contrast <5l is proportional to in 
Fourier space, and we have 



8 h (k) = M(k)0(k). 



(33) 



Throughout this paper, we do not explicitly write an argument 
of time or redshift z for simplicity in some time-dependent 
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variables. In the above equation, 5^ oc D(z) and M oc D(z) 
where D(z) is the linear growth factor. However, the potential 
<t> is the primordial one and does not depend on z by definition. 
Do not confuse with the physical Newtonian potential which 
does depend on z- The proportional factor M{k) is determined 
by the transfer function T(k) and the Poisson equation as 



2 D(z) k 2 T(k) 
M(k) = - 



3(l+z,)D(z»)H 2 Q l 



(34) 



m() 



where z* is an arbitrary redshift at the matter-dominated 
epoch. The factor (1 + z*)D(z*) does not depend on the choice 
of z* as long as z* is deep in the matter-dominated epoch, since 
D(z*) oc 1 /( 1 +z,) in that epoch. Some authors employ the nor- 
malization of the growth factor as (1 + z»)D(z*) = 1, in which 
case Eq. fl34"l i has the simplest form. 

The linear power spectrum and bispectmm of 6^ are given 

by 

P L (k) = M 2 (k)P (k) (35) 
B h (k u k 2 ,k 3 ) = M(k l )M(k 2 )M(k,)B (kuk 2 ,k 3 ), (36) 

where P<p{k) and B<p{k\ , k 2 , £3) are the primordial power spec- 
trum and bispectmm of the potential, respectively. In pop- 
ular models of non-Gaussianity, the primordial bispectra are 
uniquely related to the shape of primordial power spectrum. 

There are four models of primordial bispectra which are 
frequently considered as typical examples. Defining 



Pi = P<p(kd, 

the four models are given by the following equations: 
• the local model |42 

B^(k u k 2 , h) = 2/ NL [P X P 2 + eye] 



(37) 



(38) 



• the equilateral model [45]: 

Bf(k u k 2 , k 3 ) = 6/ NL [- (P X P 2 + eye.) - 2( J P 1 P 2J P 3 ) 2/3 

+ (P l 1/3 P 2 2/3 P 3 + 5 perm.)] . (39) 

• the folded model [46]: 

Bf (kuk 2 ,h) = 6/nl [(P1P2 + eye.) + 3{P X P 2 P^ 

- (P l 1/3 P 2 2/3 P 3 + 5 perm.)] . (40) 

• the orthogonal model j47il : 

Bf{k u k 2 ,h) = 6/nl [-3(PiP 2 + eye.) - %{P x P 2 P 3 f 13 

+3(Pi l/3 P 2 2/3 P 3 + 5 perm.)] . (41) 

With these typical models of primordial bispectrum, the lin- 
ear bispectmm Bl(^i , k 2 , £3) is given by the linear power spec- 
trum Pl(£) and a function M{k) through Eqs. ( f35T > and ( 13 61 . 
Once a model of the primordial non-Gaussianity is given, the 
non-Gaussian part of the power spectrum APx(k), AP m x(k), 
AP m {k) of Eqs. ( fTST ), ( 1301 ), ( |24"1 i and the scale-dependent bias 
of Eq. (|25l l are straightforwardly evaluated for a given model 
of bias functions. 



B. Large-scale limit of scale-dependent bias 

In the large-scale limit k — > of Eq. ( [17) . we have 
Qo(k,k') ^ Q\{k,k') * in Eqs. (Qj]i and ([19]). Accord- 
ingly, we have Qo(k) as Q\(k) as 0, and the scale-dependent 
bias of Eq. ( f25l l approximately reduces to 



Ab(k) 



2P L (kY 



(42) 



As long as the second-order renormalized bias function 
c\(k\, k 2 ) is a smooth function, we can adopt an approxima- 
tion 

c\{k',k-k')^c\(k',-k') (43) 

in the large-scale limit for the integrand of Q 2 (k). Because of 
rotational symmetry, c\{k' , - 
tude k' = \k'\, and we define 



rotational symmetry, c^{k' , -k') only depends on the magni- 



c\(k') = c\(k', -k'). (44) 

Assuming Eq. d43l for the renormalized bias function, the 
scale-dependent bias on large scales is asymptotically given 
by 



Ab(k) 



- — f ^-r&(k')B h (k,k',\k-k'\). (45) 
2P h (k) J (27r) 3 u y ' 



In this case, the scale dependence of Ab(k) is determined by 
the functional form of the primordial bispectrum B^(k, k',\k- 
k'\) in the squeezed limit k <K k' . 

For specific models of primordial non-Gaussianity intro- 
duced in the previous subsection, squeezed limits of the bis- 
pectra can be analytically derived. In the following, we as- 
sume the power-law primordial spectrum, P&(k) oc k n ^ A , 
where n s as 1 is the scalar spectral index. We define 



1 - n s 



a s = 



(46) 



to represent the deviation from the scale-free Harrison- 
Zel'dovich spectrum, « s = 1. Taking the squeezed limit 
k <sc k! in Eqs. d38l)-(l4Tli. and keeping leading orders in k/k', 
they are asymptotically given by 



Bf (k, k', \k -k'\)as Af NL Pdk)Pdk'), 

-J I -J -(l+a s ) 2 ^ 



(47) 



X P (k)P<p(k'), 

k 



Bf {k,k', \k - k'\) as 6fy L -P (k)P (k'), 
Bf (k, k\ \k - k'\) as -l2f NL tp (k)P0(k'), 



(48) 
(49) 

(50) 



where yu = k ■ k' jkk' is the direction cosine between wavevec- 
tors k and k' . Except for the equilateral non-Gaussianity, the 
lowest-order bispectra divided by Pf{k)P^(k') in the squeezed 
limit do not have explicit dependences on a s . 
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Substituting Eqs. j47l)-(l50l) into Eqs. (O and (05), we have IV. SHAPES OF RENORMALIZED BIAS FUNCTIONS 



Ab loc (k) 



2/r 



NL 



M(k) 



M(k) 



)Pdk'), 
c^(k') 



l+a s ) 



Pdk') 



Ab M (k) 



Ab° n (k) 



M(k) 
M(k) 



c^(k') 



3 J (2n 

C d 3 k' 
J (2^F 

J (2* 



) 3 k' 2 



Pdk'), 



)3 k> 



Pdk'), 



(51) 

(52) 
(53) 
(54) 



respectively. To derive the above equations, we use an ap- 
proximation Ai(\k- k'\) m Ai(k') in the squeezed limit. In the 
case of equilateral non-Gaussianity, the asymptote of Eq. (|52| i 
is somehow complicated for a general value of the spectral in- 
dex n s . In the case of scale-free spectrum, n s = 1 (a s = 0), 
Eq. ( |52| i reduces to a simpler form, 

The integrals in Eqs. d5Tli-(l55ll do not depend on k. Conse- 
quently, only the amplitudes of the scale-dependent bias Ab(k) 
are sensitive to details of biasing in the large-scale limit. The 
scaling indices are sensitive to only the primordial bispectra, 
and are independent on details of biasing. On sufficiently large 
scales, where T(k) » 1 and M(k) cc k~ 2 , the known scalings 



Ab Xoc - cc r 2 , 

Ab eql cc k°, 
Ab M cck-\ 

Ab ott - cc r 1 



(56) 
(57) 
(58) 
(59) 



hold irrespective to bias models. The scale-independence in 
the case of equilateral model holds only for the scale-free 
power spectrum, n s = 1. The above scalings are derived in 
the literature for individual models of biasing, such as the 
halo model with peak-background split [15], the high-peak 
model [19], and the local bias model l22ll . Our general argu- 
ment here shows that those scalings are general consequences 
of squeezed limits of primordial bispectra, and are indepen- 
dent on bias models. Only the proportional coefficients of 
Eqs. (T56b — (T59b depend on bias models. This finding explains 
why different models of bias in the literature give different am- 
plitudes and the same spectral index of scale-dependent bias. 

For a given model of primordial bispectrum, obtaining the 
value of scaling index of the scale-dependent bias in the large- 
scale limit is straightforward even if the bias model is not 
specified: the scaling index in that limit is just given by the 
squeezed limit of the primordial bispectrum through Eq. (|45l >, 
i.e., the large-scale behavior of Ab(k) as a function of k is de- 
termined by a combination B^ik, k ', \k-k'\)l Piik) in the limit 
of k <k k' . In the following sections, we look into the ampli- 
tudes which are determined by concrete forms of the renor- 
malized bias function of second order, c^(k\, k%). 



Detailed amplitude and shape of scale-dependent bias de- 
pend on the form of renormalized bias functions c\. In the 
lowest-order approximation of this paper, the first two func- 
tions, c\(k) and cx(ki,k2), are of primary interest. In this 
section, we derive the general form of renormalized bias func- 
tions c\i in the halo model of bias. 



A. Asymptotes of the halo bias functions 

The scale-dependent bias is previously derived in the halo 
approach with the method of peak-background split dH |H (TH 
\Y7\] . With this method, it is shown that the lowest-order term 
in scale-dependent bias, which is proportional to /nl, is pro- 
portional to the first-order Lagrangian bias parameter b 1 ^. On 
the other hand, Eulerian local bias models predict the corre- 
sponding term is proportional to the second-order bias param- 
eter £>2 |22j] . Both predictions agree with each other in the 
high-peak limit, although the difference should be important 
in physically realistic situations. Numerical simulations indi- 
cate that the corresponding term is proportional to b\, which 
agrees with the results of the peak-background split lll6l - ll8n . 

In our asymptotic prediction of Eq. (l45t . the corresponding 
term is related to the second-order renormalized bias func- 
tion, ■ Does that mean our prediction contradicts with the 
prediction of the peak-background split in the halo model? 
As we shall see below, the answer is no. To see the rela- 
tion between our general results and the prediction of peak- 
background split, one needs to derive the scale-dependence of 
the function c^(k\, kj) in the model of halo bias. 

First we introduce our notations of the halo approach. In the 
halo model, the mass M of halo is related to the Lagrangian 
radius R of a spherical cell by 



M = -n- pii R\ 



(60) 



where po is the mean matter density at the present time, or 



R = 



M 



1.162 x IQ 12 h- l M Q Q mQ 



1/3 



(61) 



where M = 1.989 x 10 30 kg is the mass of the sun. In the 
following, the radius R is always a function of a mass scale 
M through the above equation. The density variance of mass 
scale M is given by 



2 _ f d3 k 



W 2 (kR)P L (k), 



(62) 



where the function W(kR) is usually chosen to be a top-hat 
window function, 



W(x) 



3j\(x) 3 sin x - 3x cos x 



x x J 
and j\(x) is the first-order spherical Bessel function. 



(63) 
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In Ref. Ht], the large-scale asymptotes of the bias functions 
are derived in the halo approach with a universal mass func- 
tion. The mass function is said to be universal when the mass 
function has a form 



n(M)dM = §/mf(v)— , 
M v 



(64) 



where v = 6 c /cr M is a function of mass M, and S c is the crit- 
ical overdensity for spherical collapse. In the Einstein-de Sit- 
ter model, the critical overdensity is exactly independent of 
redshift, d c = 3(3tt/2) 2/3 /5 - 1.686, while it only weakly 
depends on cosmological parameters and redshift in general 
cosmology. The multiplicity function /mf(v) has a normaliza- 
tion 



JmfO) — 

Jo y 



l 



(65) 



to ensure all the mass in the universe is to be contained in 
halos. 

In the literature, several forms of the multiplicity function 
/mf are proposed. In the original Press-Schechter (PS) theory 
Hi, it has a form, 



/ps(v) 



2 
— ve 

7T 



-v 2 /2 



(66) 



Sheth and Tormen (ST) j49ll give a better fit to numerical sim- 
ulations of cold-dark-matter type cosmologies with Gaussian 
initial conditions, 



/ ST (v) = A(p) yj- 
n 



1 + 



(qv 2 )P 



\fqve 



V/2 



(67) 



where p = 0.3, q = 0.707 are numerically fitted parameters, 
and A(p) — [1 + 7r _1 ^ 2 2 _ T(l/2 - p)Y l is the normalization 
factor. When p — 0, q — 1, the ST mass function reduces to 
the PS mass function. 

Several other fitting formulas for /mf have been proposed 
with numerically improved calibrations I150l - l52ll . In Warren 
et al. lf5ltl . for example, the multiplicity function is fitted as a 
function of <x = cr M instead of v as 



/w(cr) = A (cr " + b) exp 



(-£)■ 



(68) 



where A, a, b, c are fitting parameters. The same functional 
form is applied to MICE simulations in Ref. l52ll . allowing 
the parameters redshift-dependent. Their values are given by 
A(z) = 0.58(1 +z)-° 13 , a(z) = 1.37(1 + z )-°- 15 , b(z) = 0.3(1 + 
z )-0.084 ; c(z ) = 1 036(1 + z)" 024 . As a function of v, Eq. 
can be re-expressed as 



/w(v) = MSJv) = A 



— I +h 

6, 



exp|-— |. (69) 

Oc / 



In the following, we refer to the above form as "MICE mass 
function" when the redshift-dependent parameters with MICE 
simulations are adopted. When the parameters A, a, b, c are 
redshift-dependent as in the case of MICE simulations, the 



multiplicity function explicitly depends on the redshift as 
/ MF (v, z). In this case, the mass function is not universal any- 
more. Throughout this paper, explicit dependences on time 
is notationally suppressed in the arguments of functions, and 
we adopt the notation /mf(v) even if this function explicitly 
depends on the redshift. The PS mass function is recovered 
when we formally substitute A = y/2/n6 c , a = 1, b = 0, 
c = 6 c 2 /2 into Eq. j69l . 

Using the notations introduced above, the long-wavelength 
asymptotes of the bias functions derived in Ref. Ul have the 
form, 

<£(*!,...,*„) *Z£ (M) ( |Ar £ | — > 0), (70) 
where b\;{M) is a scale-independent function defined by 

,b An), 



^(M) = |- — 



/mf(v) ' 



(71) 



and /j^p = d" fy^jdV 1 denotes nth derivative with respect to v. 

Specifically, when the PS mass function with Eq. d66l ) is 
applied, we have 



b\{M) = 



v l - 1 

6 C ' 



b\{M) = 



v 4 - 3v 2 



(72) 



and so forth, which are consistent with the results derived 
from the model of spherical collapse ll53l l54ll . When the ST 
mass function with Eq. (l67l is applied, we have 



b\{M) = i 

O c 



qv 1 - 1 + 



2p 



, 2 v 4 



1 + (qv 2 )P 
3qv z + 



2 i 2p(2qv 2 +2p - 1) 



1 + (qv 2 )P 



(73) 
(74) 



and so forth. When the MICE mass function with Eq. 
applied, we have 



a 



b\(M) 



(T- 1 + bcr a ) 
4c 2 2c a ( 4 c/c 



- a + 







a 4 



(T~ 



1 + bcr" 



(75) 



(76) 



and so forth, where cr = cr M = 5 c /v. The bias parameters 
b\(M) and b\{M) as functions of mass are plotted in Fig. [3] 
for three different mass functions considered above. 

If the asymptotes of Eq. d70b are naively used in Eq. d45l >. 
the scale-dependent bias is proportional to the second-order 
Lagrangian bias parameter Ab oc b\. Bias parameters are 
scale-independent in models of local bias. Therefore, the use 
of asymptotes of Eq. ( f70b restrict ourselves to a model of lo- 
cal bias. Local bias models predict that the scale-dependent 
bias Ab is proportional to the second-order bias parameter £>2 
in general fZM . One should note that when c^ is exactly con- 
stant, the integral in Eq. (IBTT l logarithmically diverges for the 
local model of primordial non-Gaussianity with a cold-dark- 
matter type power spectrum which has a small-scale asymp- 
tote Pl(^) 00 k ~ 3 for k — > 00. 
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FIG. 3: The scale-independent Lagrangian bias parameters derived 
from the Sheth-Tormen (solid lines), Press-Schechter (dashed lines) 
and MICE (dotted lines) mass functions for z = 0, 1,2, 4 (from bot- 
tom to top on large mass scales). The unity is added to each parame- 
ter to show the negative values in this logarithmic plot. 



The local bias model turned out not to be a good approxima- 
tion in modeling the scale-dependent bias of halos measured 
in numerical simulations with non-Gaussian initial conditions 
ifTil fl6i fl8ll . The scale-dependent bias is more or less propor- 
tional to the first-order bias parameter Ab oc b\ rather than the 
second-order parameter b\. The property Ab oc b\ is a gen- 
eral prediction of the halo approach with the method of peak- 
background split |Hl |H [Ot] . That means the method of peak- 
background split takes into account the nonlocal nature of bi- 
asing. The fact that the predictions of the peak-background 
split are better than those of local bias model in reproducing 
the results of numerical simulations implies that nonlocal na- 
ture of biasing is important to understand the scale-dependent 
bias from the primordial non-Gaussianity. 



B. Renormalized bias functions of halos without assuming 
peak-background split 

In the following several subsections, we extend the calcu- 
lation of the renormalized bias functions beyond the asymp- 
totic limit of Eq. (1701 . without resorting to the approximation 
of peak-background split. In the end, the known formula of 
scale-dependent bias derived from the peak-background split 
is exactly re-derived as limiting cases, which is one of the re- 
markable findings in this paper. 

In the literature, analytic formulas of halo bias are derived 
by more or less adopting a concept of peak-background split. 
This method is necessary because the halo approach is based 
on a statistical nature of extended Press-Schechter mass func- 



tion. In such an approach, the local mass function is obtained 
by averaging over small-scale fluctuations, while large-scale 
fluctuations are considered as background modulation field, 
which leads spatial fluctuations of number density of halos. 
Comparing the fluctuations of the halo number density field 
and those of mass, the halo bias is analytically derived. 

However, the biasing can be seen as a deterministic process 
at a most fundamental level, in which any statistical informa- 
tion is not required. One can think of getting a halo catalog 
in numerical simulations to understand the situation. Just one 
realization of the initial condition deterministic ally gives sub- 
sequent nonlinear evolutions and formation sites of halos. 

When only leading growing modes are considered in a per- 
turbation theory, any structure in the universe is determin- 
istically related to the linear density field. The biasing re- 
lation should not require statistical information of the field. 
In calculating the renormalized bias functions of Eq. (01, the 
number density fluctuations 5^ is a deterministic functional 
of linear density field 5l- Any statistical quantities, such 
as the short-mode power spectrum in the method of peak- 
background split, are not expected to appear at the most fun- 
damental level. 

From considerations above, we start from an unaveraged 
version of the original Press-Schechter (PS) formalism [48]. 
The linear density field smoothed over mass scale M, 



6 M (x) 



r (Pk 



e ikx W(kR)6 L (k), 



(77) 



is a fundamental element in the original PS formalism and its 
variants. When P(M, 5 C ) denotes the probability that (or the 
volume fraction where) the value 5 m exceeds a critical value 
6 C , the averaged number density of halos in the original PS 
formalism are given by 



«(M) = -^Ap(M,<5 c ), 
M dM 



(78) 



where n(M) is the differential mass function. In some litera- 
tures, n(M) is denoted by dn/dM, which notation we do not 
adopt. When the distribution function of 5m is exactly Gaus- 
sian, Eq. d78l is equivalent to Eq. (|64b with Eq. ([66}. 

For our purpose, we need to have a local number density 
of halos n(x, M), instead of spatially averaged one. The con- 
dition that a mass element at a particular point is included in 
a collapsed halo of mass greater than M is not stochastic but 
deterministic when the linear density field is given. Conse- 
quently, localized version of Eq. (F78b should be 



n(x,M) = 

M dM 



(79) 



where 0(x) is a step function, and x is a Lagrangian position. 
The statistical quantities do not appear in this relation. For a 
given realization of linear density field 5l, the mass element 
at a given point is either collapsed 5m > 5 C or uncollapsed 
5m < 5 C . Taking the spatial average of Eq. ( |79l , Eq. d78l of 
the original PS formalism follows, because 



(0(5 M (x)-6 c )) = P(M,6 c ). 



(80) 
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In a picture of excursion set approach ll55ll . the step function 
in Eq. d79l should be replaced by an operator yVi up (M), which 
is zero until the first up-crossing 6m > 6 C occurs as the mass 
scale M decreases from infinity, and it becomes unity below 
that mass scale. This operator yVi up is not just a single function 
of 6m, and analytic treatments are more complicated if not 
impossible. In this paper, we just use Eq. d79l in the following 
consideration for simplicity. 

It is also possible to consider the function is a general 
function which is not necessarily a step function. When this 
function is a step function, the PS mass function exactly fol- 
lows when the linear density field is random Gaussian and the 
model of spherical collapse is literally assumed. However, in 
reality, the mass function is different from the PS one even if 
the linear density field is random Gaussian. In approaches of 
universal mass function, the step function is deformed to a dif- 
ferent function to reproduce a mass function n(M) in numeri- 
cal simulations. We allow this possibility and the function 
is not necessarily a step function in the following derivation. 

The number density contrast in Lagrangian space is given 
by 6^(x) = n(x,M)/n(M) - 1 for halos of mass M, where 
Eqs. (1781 and (|79l are assumed. Taking functional derivatives 
of Eq. d79l ), we have 



6"n(x, M) 



i(ki+--+k„)-x 



2p 



5S L (ki)---S6 L (k„) 
d 



(2n) 



3« 



M 



x 



— [0 {n \5 M - 6 c W(hR) ■ ■ ■ W(k n R)] , (81) 



is a filling factor of collapsed regions. Partial derivatives in 
Eq. fl82l ) are given by 



dP(M,6 c ) = f M F(v)dlncr M 
dM 2 dM ' 

d"P(M, 6 C ) F ( "\v) 



where 



F (n\y) = — = - 



86" 



d"F 



2cr M » 



(85) 
(86) 



jn— 1 



dv" dv 



n— 1 



/mf(v) 



7=0 



In this case, Eq. ( I82t reduces to 



(-1)" 



/mf(v) d\n cr M 



F (n \v)W(k x R) ■ ■ ■ W(k„R) 



o- M " 



(88) 

To evaluate the above expression, we have useful formulas, 



d lF ln \v) 



d\ncr M \ cr M n I <T M n 



1 



/'=0 J 

dF {n \v) 



(89) 



d In cr M 



- nF ( "\v) 



An) 
Jmf 

o-m" 



(90) 



where Eq. (f77T > is used and ( - n Hx) = d"0(x)/dx" . Fourier 
transforming the above equation with respect to x and taking 
ensemble average, the renormalized bias functions of Eq. © 
or Eq. reduce to 



One can use the scale-independent parameters b};(M) = 
(-l/cr M )"/^// MF introduced in Eq. dTTT i in the above equa- 
tions. Accordingly, Eq. d88l reduces to two equivalent expres- 
sions, 



c];(k x ,...,k n ) = 



dM 



d"P(M, 6 C ) 



36" 



W(kiR) ■ ■ ■ W(k„R) 



dP(M,6 c ) 
dM 



(82) 

where R is an explicit function of M and partial derivatives of 
the window functions with respect to M do not vanish. 



C. Renormalized bias functions in universal mass functions 

Assuming a universal mass function of Eq. (l64l i, we have a 
form 



P(M,6 c ) = ^F(v), 



(83) 



cJr(Jti, . ..,*„) = ^vWifl) • • • w{k„R) 



6 C " 

A n ^(M)cr M n d 



and 



where 



6 C " 



din <jm 



W(k x R) ■ ■ ■ W(k n R) 



o-m" 



(91) 



> k n ) = b^(M)W(kiR) ■ ■ ■ W(k„R) 
A„_,(M) d 



6 C " d In cr M 



[W(hR)---W(k n R)]. (92) 



A„(M) = Y J j6 c j b l j(M) 



7=0 



(93) 



and b^(M) = 1 for consistency. There is a recursion relation, 



where 



F(v) 



r°° /mf(v) 

Jv V 



dv. 



(84) 



A n = «A„_i + 5 C "^, 



(94) 



which also guarantees the equivalence between the two ex- 
pressions of Eqs. d9TT > and (1921 . In this paper, we need the 
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FIG. 4: The functions Ai(M) and A 2 (M) derived from the Sheth- 
Tormen (solid lines), Press-Schechter (dashed lines) and MICE (dot- 
ted lines) mass functions for z = 0, 1,2, 4 (from bottom to top on 
large mass scales). 



renormalized bias functions up to second order. Explicitly, 
we have 

Ao(M) = 1, (95) 

A 1 (M) = l+d c b\(M), (96) 
A 2 (M) = 2 + 26 c b\(M) + 5 c 2 b\{M). (97) 

The functions A \{M) and A 2 (M) are plotted in Fig.|4]for three 
kinds of mass functions considered in Sec. lIV Al 

When all the magnitudes of wave vector are small, — > 
0, the window function becomes asymptotically unity, 
W(kjR) — » 1. Accordingly, the previous result of Eq. (l70t 
is correctly recovered in this limit. Therefore, Eq. (|9T1 or 
Eq. d92| > gives the form of renormalized bias functions beyond 
the approximation of large-scale limit. Scale-dependence of 
bias functions means nonlocality of the bias in general. The 
derived renormalized bias functions are scale-dependent on 
scales of mass M, which suggests that the halo bias is non- 
local on scales of halo mass, as naturally expected. Our for- 
mula, Eq. ( |9TT i or (1921 . can be used in general applications of 
iPT with the halo bias, beyond the problem of scale-dependent 
bias in this paper. 

Up to the second order, we explicitly have 



c\(k) = b\(M)W(kR) + 



1 dW(kR) 
S c d In ctm ' 



(98) 



c\(kuk 2 ) = b\{M)W(hR)W(k 2 R) 
1 + S b\{M) d 



d In ctm 



[W(MR)W(k 2 R)] . (99) 



The Eqs. d98l l and d99l complete the elements to calculate 
our prediction for the scale-dependent bias, Eq. (|231 , in the 
case of halo bias. For convenience, we define 

In(k)= f ^I n (k,k')B L (k,k',\k-k'\), (100) 



for n = 1,2, where 



2k ■ k' 

Ji (*,*') = —^W(\k - k'\R), 
k 

T 2 (k, k') = W(k'R)W(\k - k'\R). 



(101) 
(102) 



When the bias functions are given by Eqs. d98l l and (|99l , the 
functions Qi(k) and Q 2 (k) of Eqs. £[7}, (O, $M are repre- 
sented as 



Q l ik) = b\im+ 



Q 2 (k) = b\l 2 (k) + 



1 dli(k) 
5 C d In <Tm ' 
1 + 6 c b\ dl 2 (k) 
<5, 2 d\ncr M ' 



(103) 
(104) 



D. Effects of Mass Selection 

The above calculations assume that all halos have the same 
mass, M. In reality, halos of different masses are contained 
in a given sample. When the finiteness of mass range is not 
negligible, the global and local number densities of halos are 
given by 



N(M) = j dM<p(M)n(M), 
N(x, M) = j dM<p(M)n(x, M), 



(105) 
(106) 



respectively, where <p(M) is an arbitrary selection function of 
mass. When the mass of halos in a range M\ < M < M 2 are 
evenly selected, the selection function is given by <p{M) = 1 
when M\ < M < M 2 and 4>(M) = otherwise. When the 
mass range is negligibly small, <p{M') - 6d(M' - M). 

The number density contrast in Lagrangian space is given 
by S\{x) = N(x, M) /N(M) - 1 . Following a similar procedure 
to obtain Eq. ( |82"1 ), we have 



c\{k x , . . . ,k n ) 



(-1)" dM 



r 4>(M) d 
~M~dM 



d6 c 



I 



dM 



4>(M) dP(M, 6 C ) 
~M dM 



(107) 



instead of Eq. (f8Zb . Taking account of Eqs. d78l ) and d821 l. we 
find 

I dM<p(M)n(M)c];(ki,...,k n \M) 
$(k 1 ,...,k n ) = + , (108) 



dM<p(M)n(M) 
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where c^(k\, . . . , k„; M) is the renormalized bias function for 
a fixed mass M derived in the previous subsection, Eqs. ( |88t , 
(HB or d92j. Quite naturally, Eq. (11081 shows that the renor- 
malized bias functions for mass-selected halo samples are 
obtained by averaging over mass-dependent bias functions 
weighted by selected number of halos. 



E. Cancellation of highest-order bias parameters in 
Press-Schechter mass function 

In the case of PS mass function, there is a remarkable prop- 
erty that the nth-order renormalized bias function c\ only de- 
pend on next two lower-order bias parameters b^_ { and b^ iV 
as shown below. 

The scale-independent parameter b\ of Eq. (|7TT i in the case 
of PS mass function is given by 



b\(M) = 



v"-'g B+ i(v) 
6" ' 



(109) 



where H„(v) = e v2 ^ 2 (-d/dv) n e^ v '^ 2 denotes Hermite polyno- 
mials. Using a recursion relation of the Hermite polynomials, 
H„+i(y) = vH n (y)—nH n -i(v), we observe that remarkable can- 
cellations in the series of Eq. (|93T > occur and just a single term 
survives as 

A„(M) = VH„M = 5 -^bl_ x {M), (110) 
0"aT 

where b h l {M) = cr M 2 /6 c for consistency. For n = 1,2, we 
have 



A AM) 



A 2 (M) 



v 2 6 c b\(M). 



(Ill) 



Interestingly, the highest-order parameter b\ vanishes in A„, 
and it is simply proportional to the next lower-order parameter 
b\_ x in this case of PS mass function. Accordingly, the renor- 
malized bias function c\ of nth order depends only on lower- 
order bias parameters bh_^ and b\ 2 in the form of Eq. ( |9TT l 
with Eq. ( II 10b . In the next subsection, this notable property 
turns out to be the reason why the scale-dependent bias from 
the primordial bispectrum is roughly proportional to the first- 
order bias parameter b\ instead of the second-order parameter 

b L 
u 2 . 

However, the exact cancellations do not occur in other cases 
than the PS mass function. For example, in the case of ST 
mass function with the bias parameters of Eqs. ( 1731 and ( |74] >, 
coefficients of the renormalized bias functions, Eqs. ( |96l and 
(l97| > reduce to 



A\{M) = qv 2 + 



2p 



(112) 



1 + (qv 2 )P ' 

A 2 (M)^v^(M) + 2 ^ y2 ; 2 f v +1) , (113) 

1 + (qv l )P 

and so forth. The leading term of A n is still proportional to 
the next lower-order bias parameter fc L _j, but correction terms, 
which are proportional to p, additionally appear. In the case 



of PS mass function, p — and q = 1, the above equations 
agree with Eq. ( 11 lOt of n = 1 , 2 as they should. Similarly, in 
the case of MICE mass function with the bias parameters of 
Eqs. d75b and d76l ). we have 



2c 

Ai(M) = — + 1 



l+bcr"' 



A 2 (M) = -%6 c b\(M) + ~ + 2 - 



2c 



; (ic/o 2 - a + 3) 
1 + bcr a 



(114) 



(115) 



and so forth. Again, when we substitute a — 1, b = and c 
5 2 /2, the above equations agree with Eq. (II lOl i of n — 1,2. 



F. Relation to the previous formula derived by the 
peak-background split 

We explicitly show below that predictions of the peak- 
background split are re-derived as a special case of the gen- 
eral formula we have derived. First we drop the terms go, 
Qi in Eq. i25i . which are subdominant on large scales, and 
only consider the dominant term Q2(k). In this approxi- 
mation, the scale-dependent bias is approximately given by 
Ab(k)*Q 2 (k)/2P L (k),i.e., 



Ab(k) 



i r d 3 k 

2Pjk)} (2? 



^ c\(k', k - k')B L (k, k', \k - k'\). 



(116) 

Substituting the second-order renormalized bias function c 2 
in the form of Eq. d9"TT l with n = 2 into the above equation, we 
have 



where 



I(k) = 



Ab(k) 



Iiik) 

CTM 2 PL(k) 

1 

(TM 2 PL{k) 



26 2 



A 2 J(fe)+A, 



dl(k) 



d\ncr M 



(117) 



f ^W 2 (k'R)B h (k,k',\k-k'\). (118) 
J (2tt) j 



When bias parameters of the PS mass function are used, 
the parameter A\, A 2 are given by Eq. dl 1 U . and Eq. ( 11 17| i 
reduces to 



Ab{k)*Ub\l{k) + \§^. 
2 2 dm <r M 



(119) 



This equation is exactly the same as a recent prediction of 
the peak-background split lfl6l[l7ll . In the literature, the win- 
dow function W(kR) is sometimes additionally divided into 
the function I(k) because the bias is defined with respect to 
the smoothed density field in the latter case (private commu- 
nications with V. Desjacques). Earlier predictions in the litera- 
ture |Hl,[zl[l5| are reproduced from the first term of Eq. ( 11 19) . 
Although the second term is subdominant in the high- pea k 
limit, it is significant for most relevant peak heights lfl6ill7ll . 
In the local model of non-Gaussianity, the second term van- 
ishes on sufficiently large scales. 
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Importantly, our derivation does not use the approxima- 
tion of the peak-background split to obtain the same pre- 
diction. Therefore, the known result of Eq. (II 191 turns out 
not to particularly depend on the approximation of the peak- 
background split. However, we have found that this equation 
is consistent only with the PS mass function. That is not sur- 
prising, because the arguments of peak-background split use 
an auxiliary Gaussian field as an fundamental field in the anal- 
ysis. The original PS mass function is properly derived by 
assuming the Gaussian statistics for the linear density field. 
Applying the Gaussian statistics to the auxiliary field is only 
consistent with the original PS mass function in the peak- 
background split. 

Our derivation, on the other hand, does not use any aux- 
iliary Gaussian field, and therefore gives natural extensions 
to the formula with general mass functions as we have seen. 
Even when subdominant terms on large scales, Qo and Q\, 
are neglected, our result predicts that the known result of 
Eq. ( 1119b should be modified to Eq. (1117b with Eqs. d9oT l and 
(13, i.e., 



Ab(k) * ^ [(2 + 26 c b\ + 5 2 b\) I(k) 



(l+6 b\) 



dl(k) 



(120) 



Taking the ST mass function for example, the above equation 
reduces to 



Ab(k) 



qS c b\ 1 p(qv 2 + 2/7 + 1) 



1 + (qv 2 )P 
1 P 



I(k) 



1 + (qv 2 )P 



dl(k) 
din (Jm 



(121) 



This equation is still a new result in the literature, even though 
other correction terms, Qo and Q\, are neglected. The previ- 
ously known result, Eq. dl 19b , is obtained only when the PS 
mass function is assumed, p — 0, q — 1. Taking the MICE 
mass function, the coefficients of Eq. ( 1121b are replaced by 
Eqs. (Q3I) and ( fTT3T l instead of Eqs. (fTTSl and ( fTBT l. 

One should note that although the ST (MICE) mass func- 
tion is derived from Gaussian simulations, contributions of the 
primordial non-Gaussianity to the mass function are higher or- 
ders in Eq. (1121b . and can be neglected in the lowest-order ap- 
proximation of this paper. Although higher-order corrections 
are beyond the scope of this paper, they could be relevant in 
the quantitative analysis of actual data. 

In the specific models of primordial non-Gaussianity, large- 
scale limits of the function I(k) with asymptotic bispectra of 



Eqs. (I47li-(l50b are given by 
4/nl 



I loc (k) 



M(k)' 



k "Jl+la, Z 72 



M(k) 



where 



y a (M) 



1 r d 3 k 
cr M 2 J (2n) 



■k- a W z (kR)P L (k). 



(122) 
(123) 
(124) 
(125) 

(126) 



Predictions of the scale-dependent bias in these specific mod- 
els are made when the above equations are substituted in 
Eq. (1121b . In the scale-free power spectrum, n s = 1 {a s = 0), 
Eq. (1123b is simply given by 



r ql (k) 



8/f 



NL ,2 



M(k) 



k £ y 2 , (n, = 1). 



(127) 



For highly biased objects, A 2 ~ qv 2 5 c b\ » Au the domi- 
nant term of Eq. (1121b is given by 



Ab(k) 



qd c b\ 



I(k), 



(128) 



which agrees with previous predictions ]5j |6j [15|] in the case 
of PS mass function, q — 1 . In the case of ST mass func- 
tion, we find the factor q 0.707 should additionally be 
present. Since the ST mass function gives better fit to the halo 
mass function in numerical simulations, our result suggests 
that these previous predictions overestimate the amplitude of 
scale-dependent bias in the high-peak limit. In the case of 
MICE mass function, the coefficient q in Eq. (1 1 28b is replaced 
by 2c 15 2 = 0.729(1 + zf QmA . 

In fact, recent numerical simulations actually prefer that 
the previous theoretical predictions overestimate the ampli- 
tude of scale-dependent bias B23I427I1 . When the parameter q 
in Eq. ( 1128b is freely fit to numerical simulations, they found 
q — 0.6 - 1 although the value varies from simulation to sim- 
ulation, depending on the algorithm to identify the halos, etc. 
However, it is encouraging that the value is not so different 
from q = 0.707. It is also natural that numerical simulations 
do not fit well to the high-peak formula of Eq. ( 1128b , because 
the correction terms in Eq. (1121b are not negligible for rele- 
vant ranges of halo mass. More quantitative comparisons of 
the newly derived Eq. ( 1120b with numerical simulations, to- 
gether with estimating higher-order corrections and improv- 
ing nonlocal bias models, are left for future work. 



V. NUMERICAL COMPARISONS 

In this subsection, we numerically evaluate the analytical 
results derived above, and compare them with other approxi- 
mate methods. For this purpose, we use the halo bias as a spe- 
cific example, where the bias functions are given by Eqs. 
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FIG. 5: The functions Q n (k) divided by /nl^l(^) (dashed: n = 0, 
dotted: n = 1, solid: n = 2) at redshift z = 1, Negative values 
are shown in thin lines. The halo model with mass M - 10 14 M o is 
assumed in calculating Q\{k) and Qzik). Different panels correspond 
to different models of primordial non-Gaussianity as indicated. 



FIG. 6: The functions Q„(k) divided by /nl^l(^) with a fixed scale 
k = 0.005 h Mpc -1 at redshift z = 1, plotted against the halo mass 
(dashed: n = 0, dotted: n = 1, solid: n = 2). Negative values are 
shown in thin lines. Different panels correspond to different models 
of primordial non-Gaussianity as indicated. 



and ( l99l > up to second order. We assume the halo sample has 
a redshift of z = 1, as a typical example. 

For numerical integrations, it is useful to rewrite the func- 
tion Qo of Eqs. ( fPTT i and (11 St . and the functions I\{k) and 
of Eqs. (1 1 00] > — (1 1 02b into the following two-dimensional 
integrals: 



Qo(k) = 



4k 2 



XOO y-»l 
dr I dx 



2rx — 



4 r 2 (l -x 2 ) 



7 1 + r 2 — 2rx 
xB L (k,kr,ky/l +r 2 -2rx), (129) 



I 1 (k)=—J dr J dxrxW{kR^\ +r 2 -2rx) 



Ii(k) = 



xB L (k,kr,ky/l + r 2 -2rx), (130) 
— J dr J dx r 2 W(krR)W (kR^ll+r 2 - 2rx) 
xB L [k,kr,k^l +r 2 -2rx). (131) 



In terms of the functions I\{k) and the functions Q\{k) 

and Qiik) are given by Eqs. (1 1 031 > and ( 1104b . In numerically 
calculating the derivatives dI n /dlncr M , it is useful to note 
that the derivative of the top-hat window function [Eq. ( |63l l 
is given by 



dW(x) 3 j2{x) 3(jc 2 - 3) sin x + 9x cos x 
dx x x 4 



(132) 



where jz(jx) is the second-order spherical Bessel function. 

In Fig. [5] the functions Q„{k) (n = 0, 1,2) are plotted in 
the case of z = 1 and M = 10 14 M Q . We assume four models 
of primordial non-Gaussianity, Eqs. (|38)-(|4T), and bias func- 
tions with ST mass function, Eqs. {73), (73), ® and d99l ). 
As we have seen in Sec. IIII Bl the function Qi(k) is dominant 
over Qo(k) and QiW on sufficiently large scales. In Fig. [6] 
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FIG. 7: Scale-dependent bias Ab{k) divided by /ml for halos of mass 
M = 10 14 M o . Solid lines correspond to predictions of iPT with ST 
mass function. Other lines correspond to predictions of the peak- 
background split (PBS, dotted lines), the low-A: limit of iPT with ST 
mass function (dashed lines), and the high-peak limit of PBS (dot- 
dashed lines). Positive and negative values are represented by thick 
and thin lines, respectively. When we use the PS mass function for 
the iPT, low-fc limit of iPT and PBS gives the same results. 



the same functions are plotted against the halo mass with a 
fixed scale k = 0.005 h Mpc -1 . On large scales, the function 
Qi{k) is dominant as expected, and the scale-dependent bias 
is approximately given by Eq. (l42l . 

In Fig. |2l the scale-dependent bias Ab{k) divided by /nl 
is plotted for M = 1O 14 M with four models of primordial 
non-Gaussianity. In solid lines, the predictions of iPT with 
ST mass function [Eqs. (25]), <|98]), <|99]), (73]) and O] are 
plotted. The predictions of the peak-background split (PBS) 
[Eq. dl 19H are plotted in dotted lines. In dashed lines, the 
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FIG. 8: The same as Fig. [7] but the scale-dependent bias Ab(k) is 
normalized by the large-scale, high-peak approximation Ab Hr (k). 



predictions of the large-scale (low-fc) limit of iPT with ST 
mass function [Eq. ( I121H are plotted. In dot-dashed lines, 
predictions of the PBS with high-peak and large-scale limits 
are plotted. The last predictions are given by 



Ab HP (k) 



M(k) 



(2, (local), 

4k 2 j2 (equilateral), 

3ky\ (folded), 

—6kj\ (orthogonal), 



(133) 



which are derived from Eqs. ([122), (TT271 (fT24l , ([T251 and 
the first term of Eq. ( 11 19) . For the equilateral model, we use 
the approximation n s — 1 for simplicity, and adopt Eq. ( 1127t 
instead of Eq. ( 1123l l. The first case of Eq. (11 33b is exactly 
the same as the original formula for the local model derived 
in Refs. JUI^]. The third and forth cases exactly match the 
results of Ref. lfl5ll for the folded and orthogonal models. For 
the equilateral model, however, the prefactor4 is replaced by 6 
in Ref. [15], because angular-dependences of the bispectrum 
are neglected in the latter. The factor 4 is more accurate as 
shown below. In Fig. [8j relative amplitudes with respect to 
Ab HP (k) of Eq. (TT33T > are plotted. 

As expected, the low-k limit of iPT (using only Q2) 
is an extremely good approximation on large scales (k < 
0.05 n/Mpc). The large-scale asymptotes of iPT with ST mass 
function is smaller than those of peak-background split. These 
differences originate from the different mass functions. If we 
adopt PS mass function in iPT, the asymptotes of iPT and PBS 
exactly agree. 

The simple predictions of Ab HP (k) have correct slopes in 
large-scale limits. The amplitudes are not so accurate for 
non-local-type non-Gaussianities, even if we assume PS mass 
function instead of ST mass function. On smaller scales 
(k > 0.01 n/Mpc), the slopes of Ab np (k) are not accurate 
enough. The deviations from the peak-background split in 
k > O.Oln^'Mpc are also pointed out in earlier work in the 
high-peak limit of thresholded regions §Q [56|]. However, 
higher-order loop corrections which we do not consider in this 
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FIG. 9: The normalized scale-dependent bias in large-scale limits 
(k — > 0, however, we numerically adopt k = 10~ 4 /i/Mpc to plot 
this figure). As in Fig. [8] the values are normalized by the simple 
approximation Ab HP (k). Meanings of different types of lines are the 
same as in Figs.[7Jand[8] 



paper might affect the behaviors of scale-dependent bias on 
scales of k > 0.05/i _1 Mpc. 

In Fig. [9] the normalized scale-dependent bias against the 
halo mass is plotted in the large-scale limits (k — > 0). The 
purpose of this figure is to illustrate the mass dependence of 
the plateaux on large scales in Fig. [8] The predictions of peak- 
background split and high-peak limit are asymptotically agree 
with each other for high-mass halos, as they should. In large- 
mass limits, the predictions of peak-background split and iPT 
are different by a factor of q ^ 0.7, which is explained by 
Eq. (EE)- 



VI. SCALE-DEPENDENT BIAS IN REDSHIFT SPACE 

Generalizing the results obtained so far in real space to 
those in redshift space is fairly straightforward in the frame- 
work of iPT. The redshift-space distortions of scale-dependent 
bias in the presence of local-type non-Gaussianity have been 
investigated in the literature II57| - |59|| . relying on the high-peak 
limit or the peak-background split. With the iPT, it is straight- 
forward to generalize those results without restricting to the 
local-type non-Gaussianity, and we do not need to rely on the 
high-peak limit or the peak-background split. 

When the weak time dependences in the perturbation ker- 
nels are neglected, the nth order perturbations are approxi- 
mately proportional to D", and the redshift-space counterpart 
of the kernel L„ is given by rf36tl 

K = R (n) L n , (134) 

where /? (n) is a 3 x 3 matrix with elements 



(135) 



where / = din D/d In a = D/HD is the linear growth rate, 
and Zi is the z'th element of the unit vector z along the line of 
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sight. The perturbation kernels L„ in the multipoint propaga- 
tors of Eqs. (|5]l and (0 are replaced by U n in redshift space. 
Specifically we have 



Eventually, the non-Gaussian part of the power spectrum in 
redshift space is given by 



k ■ L s n = (k + nffxkz) ■ L n , 



(136) 



where p = k ■ z/k is the direction cosine of the wavevector k 
with respect to the line of sight. 

Substituting the resulting multipoint propagators in red- 
shift space into Eq. (fT3l . we have an expression for the non- 
Gaussian part of the power spectrum in redshift space. In eval- 
uating the expression, the following integrals are useful: 



/ 
/ 

/ 
/ 



^L 2 (k', k - k')B L (k, k',\k- k'\) = 1 ±R 2 (k), 
(2ny 7 k L 



(137) 



^-L u {k')Uj{k - k')B L (k, k', \k - k'\) 
(Zny 



k 4 



2Ri(t) - -R 2 (k) 



~^Mk), (138) 



^-^(k'^ik - k')B h (k,k', \k - k'\) = ~Qi(k), 
(27T) 1 2 k 



d 3 k' 
(2n) 



(139) 

: c\(k', k - k')B L (k, k', \k - k'\) = Q 2 (k), (140) 



where the functions Q\{k) and Q 2 (k) are defined by Eqs. ( fTTT i, 
( fT8l . ( TT9b . and the functions R\{k) and R 2 (k) are defined by 



Rn(k) 



r d 3 k' . ,. 

TR n (k,k' 

J (2nf 



)B L (k,k',\k-k'\), 



where 



Ri 

R 2 



kk' 

k 2 k' 2 - (k ■ k'f 
k' 2 \k-k'\ 2 ' 



(141) 

(142) 
(143) 



Eqs. ( 1 1 37b — ( fT40b are derived by noting that integrals on the 
left-hand sides are functions of only k, and using the rota- 
tional covariance (similar technique is used in Ref. [36i;]). The 
function Qo(k) is related by 



Q (k) = 2R 1 (k)--R 2 (k). 



(144) 



The functions Ri(k) and R 2 (k) reduce to two-dimensional in- 
tegrals as 



/? 1 (fe) = _J dr J dx rx B L (k, kr, k Vl + r 2 - 2rx) , 



(145) 



Ki /-ex. pi 2(1 _ 



2n _ x 2- ) 
2rx 



B L (k,kr,kil+r 2 -2rx). (146) 



AP s x (k,p) = [h(k)+fp 2 ] |2(1 +ff) 2 R l {k) 
[4.. 1 



R 2 (k) 



+ d+fp 2 )Q 1 (k) + Q 2 (k)\. (147) 



When we put / = in the above equation, the power spec- 
trum in real space, Eq. (fT~5T >. is recovered, In the large-scale 
limit k -> 0, we have Ri(k),R 2 (k),Qi(k) <sc Q 2 (k), and the 
asymptotic form of the power spectrum is given by 



AP x {k,p) * [bm + fp 2 ] Qi(k). (148) 



The angular dependence of the power spectrum is conve- 
niently decomposed into multipoles as 



APyiKfx) = Y J Ap l (k)P,(jx), 



(149) 



where Piip) are the Legendre polynomials, and Api(k) are the 
multipole moments. Using the orthogonal relations of the 
Legendre polynomials, we have 



Ap,(k) = 



21+ 1 



i: 



dpPi(p)AP x (k,p). 



(150) 



The monopole component Apo(k) corresponds to the angular 
average of the power spectrum in redshift space, and higher- 
order moments characterize anisotropics in the power spec- 
trum, relative to the line of sight. We decompose each mul- 
tipole moment as pi(k) = pf(k) + Api(k), where pf(k) is the 
multipole moment of the Gaussian part of the power spectrum 
Pf(k,fi). 

The multipole decomposition of the linear power spectrum 
in redshift space P G (k,p) = [b\{k) + fp 2 ] 2 ?^) is well known 
MM- 



p°(k) = j[£i(£)] 2 + |fciC*)/ + l -f \ P\(k). 
p^(k) = (^b 1 (k)f+jf 2 \p L (k), 



p C A (k) = ^fPdk). 



(151) 
(152) 
(153) 



For the multipole decomposition of the non-Gaussian part, 
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Eq. ( 11501 ). we have 
Ap (k) = 2 



10" 



/ 2f 2 f I 2f f 2 , 



-2 



2/ 4/ 2 3/ 3 l 2 4/ 2Z 2 

1 1 h — I 1 b\(k) 

21 35 35 \7 21 15 1 



Ri(k) 



Apiik) = 4/ 

-/ 
+ 2/ 



G2(fe), 

5/ 2 



(154) 



1 4/ 

3 + T + lT + l3 



7 



*i(Jfc) 



8 32/ ll/ 2 /16 13/\, 
— + — - + — — + — + — - \b) 
21 49 21 \2\ 21 / 



(k) 



l - + 2 -t + hm 

3 7 3 



Q,{k)+ 1 ^Q 2 {k), 



Ri(k) 
(155) 



/f/> 4 (&) = 16/ 2 



2 3f 1 

35 + 77 + 35 W 



*i(*) 



1Z! 

35 

32/ 3 



16 26/ 

Y + n +b](k) 



R2(k)+^-Qm, (156) 



(157) 



Although the above expressions are messy, most of the terms 
are negligibly small in the large-scale limit k — > 0, and only 
terms proportional to Q 2 (k) are dominant, i.e., 



4po(k) ~\ J -+bi\ Q 2 (k) 



A P2 {k) 

Ap 4 (k),Ap 6 (k) «: Ap (k),Ap 2 (k), 



^-Q 2 (k), (158) 
(159) 



where the linear (Gaussian) bias parameter b\ is generally 
scale-independent in the large-scale limit. The factor f/3 + b\ 
in the monopole component is previously derived in a spe- 
cial case of the local-type non-Gaussianity with the high-peak 
limit 

When the bias is large enough, » c\ » 1, and Q 2 » 
Qi » R i , R 2 , we have 



[bmfp^k) » P $(k) » pZ 



(160) 



Ap Q (k) * bmQiik) » ^ 2 (fe) » » 4> 6 (*)- (161) 

The above equations show that the redshift-space clustering 
reduces to the real-space clustering when the bias is large 
enough. This fact is naturally expected because the peculiar 
velocities of high peaks or high-mass halos are sufficiently 
small. 

In Fig. [10] relative differences of the lower-order multi- 
poles, Api{k)l p^(k), are plotted for M = 1O 14 M . We only 
plot the lowest two multipoles I — 0, 2, since the higher-order 
multipoles / = 4, 6 are small enough. The scale dependences 
of multipoles with I = 0, 2 have a similar slope on large scales 
but different amplitudes. These properties are explained by 
the dominant contributions on large scales, Eq. (1158) . The ra- 
tio of those two multipoles, in the large-scale limit, is given 
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FIG. 10: Relative differences of the multipoles in the halo power 
spectrum Api(k)/pf(k) with M = 10 l4 M o (dashed lines: / = 0, dotted 
lines: / = 2). Solid lines corresponds to the relative differences of the 
real-space power spectrum, which are the same as the solid lines in 
Fig-El Positive and negative values are represented by thick and thin 
lines, respectively. 



by 



A P2 (k) 2/ 



A PQ (k) f + 3bi 1 + 3/0' 



(162) 



where = f/b\ is the linear redshift-space distortion param- 
eter. This ratio is small for highly biased objects, or massive 
halos. 

Relative monopole components Apo(k)/p^(k) on large 
scales are not so different between real space and redshift- 
space. This property is true even in the low-mass halos. 
Therefore, the scale-dependent bias is not so affected by the 
redshift-space distortions when only the monopole component 
is concerned. However, the quadrupole component may be 
used in constraining the primordial non-Gaussianity. 

In Fig.QTl relative differences of the lower-order multipoles 
are plotted against the mass of halo, where the scale is fixed to 
k = 0.005 /i^'Mpc. As described above, the relative monopole 
components are similar in real and redshift spaces, irrespec- 
tive to the mass of halos. They are asymptotically the same in 
the high-mass limit, because of Eqs. ( 1160t and ( 11611 >. In the 
same limit, the ratio Ap 2 /Apo approaches to 20/3 according 
to Eq. ( fT62l . 



VII. CONCLUSIONS 

The iPT is a general framework of the perturbation theory 
in the presence of bias. In this paper, we first apply this frame- 
work to deriving the relation between the scale-dependent bias 
and the primordial non-Gaussianity. Approximations such 
as the peak-background split and the high-peak limit, which 
are usually adopted in the literature to estimate the scale- 
dependent bias, are not required. The redshift-space distor- 
tions of the scale-dependent bias are also evaluated. Thus, in 
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FIG. 11: Relative differences of the multipoles in the halo power 
spectrum as a function of the halo mass for a fixed scale k = 
0.005 /r'Mpc. Meanings of different types of lines are the same as 
in Fig. [TO] (solid lines: real space, dashed lines: / = 0, dotted lines: 
1 = 2). 



this paper, we have derived the most general formula so far of 
the scale-dependent bias with primordial non-Gaussianity in 
the literature. 

For the scale-dependent bias in real space, the most funda- 
mental equation in this paper is provided by Eq. ((25}, where 
Q n (k) is linearly dependent on the primordial bispectrum. We 
find that the slope of the scale-dependent bias in the large- 
scale limit is determined only by primordial bispectra in the 
squeezed limit, and is independent on detailed models of bias. 
This property explains the fact that different models of bias 
have predicted the same slope of the scale-dependent bias in 
the literature. 

We derive the shape of renormalized bias functions, gener- 
alizing the concept of simple Press-Schechter approach. The 
general expression of renormalized bias functions in this ap- 
proach is given by Eq. (l82l . In the case of universal mass func- 
tion, the renormalized bias functions are given by Eq. (|9H , or 
equivalently Eq. j92l . The previously known results in the 
approximation of peak-background split are reproduced when 
the PS mass function is assumed in our results [Eq. ( 11 19H . 
When the mass function deviates from the PS form, our results 
suggest that the form of scale-dependent bias should be cor- 
rected. The general formula of scale-dependent bias on large 
scales is given by Eq. dl20t . This equation is one of the main 
outcomes in this paper. Most of the previous results regarding 
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the scale-dependent bias from the primordial non-Gaussianity 
are derived as special cases of this equation. 

The evaluations of the redshift-space distortions in the 
scale-dependent bias are straightforward in the framework of 
iPT The result is given by Eq. (11471 1. or in terms of mul- 
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While highly biased objects have large amplitudes of power 
spectrum, the number of objects is small and the shot noise is 
large. Thus highly biased objects are not suitable for testing 
the primordial non-Gaussianity. On the other hand, the ampli- 
tude of power spectrum is small for less biased objects, and 
the clustering signals are small. Consequently, there should 
be an optimal objects with sufficiently large bias and suffi- 
ciently large numbers at the same time for realistically con- 
straining the primordial non-Gaussianity by galaxy surveys. 
The high-peak limit or the peak-background split are not nec- 
essarily valid in some cases. The results of this paper pro- 
vide the most accurate formula of the scale-dependent bias in 
the literature. They should be useful in theoretical investiga- 
tions as well as in constraining the primordial non-Gaussianity 
with realistic galaxy surveys. Applications of the results in 
this paper, including the Fisher analysis of the future galaxy 
surveys, higher-order analyses of primordial non-Gaussianity, 
are now in progress. For more accurate modeling of the scale- 
dependent bias, it should be necessary to improve the nonlocal 
bias model beyond the simple halo approach. Investigations 
in this direction will be addressed in future work. 
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